Pairing, off-diagonal long-range order, and quantum phase transition in strongly 
attracting ultracold Bose gas mixtures in tight waveguides 



o 
o 



M. D. GirardeaiEl 

College of Optical Sciences, University of Arizona, Tucson, AZ 85721, USA 

(Dated: May 5, 2009) 

A model of two ID ideal Bose gases A and B with strong AB attractions induced by a p-wave AB 
Feshbach is studied. The model is solved exactly by a Bose-Bose duality mapping, and it is shown 
that there is no A-component or B-component Bose-Einstein condensation and no AB-pair ofT- 
diagonal long-range order (ODLRO), but both AA-pair and BB-pair ODLRO. After generalization 
by adding even-wave AA and BB repulsions and reducing the strength of the odd-wave AB attraction 
by Feshbach resonance detuning, a quantum phase transition occurs between a phase with AB 
contact nodes and one with no such nodes. 
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Strong interatomic interactions and correlations occur 
in ultracold gases confined in de Broglie waveguides with 
transverse trapping so tight that the atomic dynamics is 
essentially one-dimensional (ID) with confinement- 
induced resonances [l|, 0] allowing Feshbach resonance 
tuning 3] of the effective ID interactions to very large 
values. This has led to experimental verification [J, 0, [a| 
of the fermionization of bosonic ultracold vapors in such 
geometries predicted by the Fermi-Bose (FB) mapping 
method 7], an exact mapping of a ID gas of bosons with 
point hard core repulsions, the "Tonks- Girardeau" (TG) 
gas, to an ideal spin-aligned Fermi gas. The "fermionic 
Tonks-Girardeau" (FTG) gas [8, 9], a ID spin-aligned 
Fermi gas with very strong attractive interactions, can 
be realized by a 3D p-wave Feshbach resonance as, e.g., 
in ultracold ''"K vapor [l^. It has been pointed out 
d, 13, [^, [llj that the FB mapping 7] can be exploited 
to map the FTG gas to the ideal Bose gas. The very 
strong fermion-fermion attraction in the FTG gas leads 
to fermion pairing and superconductive off-diagonal long- 
range order (ODLRO) of the two-fermion density matrix 

m- 

There are a number of models of strongly interacting 
ID ultracold gas mixtures which are exactly soluble by 
generalizations of the FB mapping In [Hill the 

properties of one such soluble model, a mixture of a TG 
Bose gas A and an ideal Fermi gas B, with point hard- 
core AB interactions, were investigated in detail. Here 
I shall discuss another model, a mixture of two differ- 
ent ideal Bose gases A and B, with an AB interaction 
of FTG form. This model is exactly soluble, and it will 
be shown that it has very unusual behavior: The strong 
AB attraction destroys the ground state Bose-Einstein 
condensation (BEG) and single-particle off-diagonal long- 
range order (ODLRO) of both components A and B, and 
it induces both AA and BB pairing manifested in super- 
conductive ODLRO of both the two-A and two-B density 
matrices, although there are no AA or BB interactions. 
Furthermore, there is no AB pair ODLRO in spite of the 
strong AB attractions. It will also be shown that if the 



AB attraction is a finite odd-wave attraction rather than 
the infinite FTG limit, and there is also a repulsive even- 
wave AB interaction of Lieb-Liniger (LL) delta function 
form [isl , then there is a quantum phase transition as the 
coupling constants are varied, between a phase in which 
there are no AB contact nodes and only the repulsive 
LL interaction acts, and another phase in which there 
are AB contact nodes and only the attractive FTG-like 
interaction acts. 

FTG interaction and nodal structure: The FTG gas 
is a spin-aligned ID Fermi gas with infinitely strongly 
attractive zero-range odd-wave interaction induced by a 
p-wave Feshbach resonance. It is the infinite ID scat- 
tering length limit oid — > —oo of a ID Fermi gas with 
zero-range attractive interactions leading to a ID scatter- 
ing length defined in terms of the ratio of the derivative 
^I^ of the wave function to its value at contact: '^{xjk — 
0+) = -^{xjk = 0-) = -ai_D*'(xjfe = 0±) where the 
prime denotes the derivative with respect to xjk 0, S, @] ■ 
The FTG limit aiu —oo corresponds to a ID zero- 
energy odd-wave scattering resonance reachable by Fes- 
hbach resonance tuning to a ID odd-wave confinement- 
induced resonance [H, 0, Q . There are several different 
zero-range pseudopotentials which generate this contact 
condition. One representation |8|| is Vo = go6 {xjk)d± 
where d±'li{xjk) = (l/2)[*'(0-h) -f *'(0-)]. The FTG 
limit aiD — )■ — oo is equivalent to go +oo. This 
representation explicitly exhibits an odd-wave projection 
property of the interaction, i.e., it vanishes on even func- 
tions of Xjk- 

Since there is no particular symmetry under exchange 
of particles of the different species A and B, AB scat- 
tering in all partial waves is possible, but usually s-wave 
scattering dominates. However, in the neighborhood of a 
p-wave AB resonance, p-wave scattering dominates, and 
gives rise in ID to an odd- wave effective interaction of 
FTG form 0, S, i]. For this it is important to real- 
ize that the odd-wave projection property of Vo requires 
only local antisymmetry, i.e., it guarantees that if Xi is 
an A-particle position and yj a B-particle position, then 
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Vo^{xi,yj) is nonzero only if ^ has a node at Xi = yj, 
where "^{xi,yj) = — ^'(yj, a;^) in an infinitesimal neigh- 
borhood of the node. It is not necessary that there be 
global antisymmetry under such exchange, and indeed, 
the ground state derived in the following sections has 
no such global antisymmetry. In the presence of both 
a ID even- wave AB Lieb-Liniger interaction geS{xi — yj) 
flit generated by 3D s-wave AB scattering and a ID odd- 
wave resonance generated by a 3D p-wave Feshbach reso- 
nance for AB scattering, the wave function can lower its 
energy by developing nodes at Xi = yj so as to kill the 
repulsive even-wave interaction and activate the strong 
ID odd-wave attraction. 

The contact discontinuities of 5* [11] can also be un- 
derstood as a zero-range limit xq 0+ and Vq — > oo 
of the two-body scattering solution for a square well of 
width 2xo and depth Vq, where the limit is carried out 
such that VqXq approaches a finite, nonzero limit @, Q. 
In the untrapped case, the exterior solution for scatter- 
ing length aAB —> —oo is constant (-1-1 for x — y > 
and -1 for x — y < 0), and the interior solution is 
sin[K(x — y)] with k ~ ^j2^iVa/h? = tt/2xq and fi the 
effective mass mAmBimA + ms)', hence kxq — 7r/2 and 
VqXq = {wh)'^ /Sfi. In the zero-range limit the interior 
kinetic energy — s- +oo and potential energy — > — oo, but 
their sum remains zero, the ground state energy. Since 
the FTG interaction acts only on odd waves, the FTG 
interaction vq should be written as vq — vPq where v is 
the above square well and Pq is the odd-wave projector. 

Many-body ground state: Assuming trapping on a ring, 
the Hamiltonian consists only of the kinetic energy oper- 
ators of components A and B plus the AB FTG interac- 
tion: 
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where Vo is the previously defined odd- wave FTG interac- 
tion. The wave functions satisfy periodic boundary con- 
ditions with periodicity length L (the ring circumference) 
with respect to all the Xi and yi. The scattering length 
uab is — oo in the FTG limit. For square well width 2a;o 
nonzero but very small, the unnormalized ground state 
^'o is constant (say ±1) when all \xi — yj\ > 2x{)^ ex- 
cept for sign changes as each Xi — yj varies from ~xq 
to Xq in accordance with the internal wave function 
± siii[K,{xi — yj)]- The condition k, = yj2^VQ/h? = tt/2xq 
determines the well depth Vq such that the scattering 
length OAB is — oo, and the FTG limit is obtained by 
letting — > and Vq ^ oo in accordance with this 



condition. In this limit the internal wave function be- 
comes invisible and V&o appears to jump discontinuously 
between ±1 whenever an A-particle passes a B-particle, 
but there are hidden nodes at Xi — y j = at the centers 
of the wells. The ground state energy Eq=^ O^eneraliz- 
ing the situation for the pure FTG gas 0, H, II [H . *o 
maps to a "model state" 5' mo consisting of two noninter- 
acting ideal Bose gases totally Bose-Einstein condensed 
into their ground orbital, which is a trivial constant for 
periodic boundary conditions: '^q = ^'mo-^^ = M and 
MO = 1 where M is the mapping function 

Na Nb 

M{xi, ■ ■ ■ ,XNA;yi, ■ ■ ■ .VNb) = WWs&^{xi - yj) (2) 

1=1 ]=i 

where the sign function sgn(a;) is 4-1 (—1) if a:; > (a; < 
0). Although the mapped bosonic state ^ mo is a trivial 
constant outside the square wells, the interior wave func- 
tion ± sin(K|a;i — yj\) vanishes with cusps at Xi — yj = 0. 
Therefore, physical consistency requires the presence of a 
zero-diameter hard core interaction added to the square 
well. The mapped Bose gas is then not truly ideal, but 
rather a TG gas with superimposed attractive well, whose 
nontrivial interior wave function becomes invisible in the 
zero-range limit, simulating a mixture of two noninter- 
acting ideal Bose gases insofar as the energy and exterior 
wave function are concerned. The densities of compo- 
nents A and B are trivial constants in the ground state, 
but the off-diagonal elements of the reduced density ma- 
trices of ^0 are quite nontrivial and interesting, due to 
the effects of the discontinuities in M; they will be dis- 
cussed in later sections. 

The periodic boundary conditions impose constraints 
on the values of Na and Nb- Suppose that the positions 
of all B-bosons and all but one A-boson, say Xi , are fixed 
on the open interval (0,L). Starting with that A-boson 
at xi > but to the left of all the other particles and 
moving it to a position xi < L but to the right of all 
the others, one counts Nb sign changes of due to the 
AB contact nodes, so Nb must be even for L-periodicity 
in the A-boson coordinates. Repeating this process with 
all A-bosons and all but one B-boson fixed and moving 
that B-boson instead, one counts Na sign changes; hence 
Na must be even for L-periodicity in the B-boson coor- 
dinates. 

One-particle density matrices and momentum distri- 
butions: Generalizing the derivation in [l^ . [l^ [l7| , one 
finds that the one-particle density matrix of component 
A is 



J 



piAix,x') = NaL 



-Na-Nb 



'^oix,X2,--- ,xna;Y)^o{x\x2, - ■ ■ ,XNA;Y)dx2 ■ ■ ■ dxNAdY ^ nA[I{x,x')]^'' (3) 
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where Y = (yi,-" iJ/Wb), "-a = Na/L, and I{x,x') = 
-^~V-L/2Sgn(a;-C)sgn(a;'-0rfC = l-2|x-a;'|/L. In the 
thermodynamic Hmit where Nb — > oo and L —^ oo such 
that Ng/L ng with finite and nonzero number den- 
sity np ng, one has [/(x, x')]^^ ^-2nB\x-x \ anal- 



ogy with the one-component FTG gas case [l2, 1^ 12 1- 
Hence pia{x,x') ~ n^e^^"^'^^^ I, and by interchange 
of A and B piB{y,y') = nse"^"'*''^"^ L Their Fourier 
transforms UkA and UkB, normahzed to J2k''^kA = Na 
and rifes — Ng where the allowed momenta are 
k ~ viiT / L with V = 0, ±1, ±2, • • • ), are Lorentzian 
discrete momentum distributions UhA — ."^"/"f, and 

nkB = "fc2 ■ The Fermi-like Lorentzian shapes are 
strong modifications of the ideal Bose gas distributions 
NASkO and NgSkO due to the infinite AB attraction; no 
trace of BEG of components A and B, and of the asso- 
ciated ODLRO of piA and pig, remains. Nevertheless, 
at fc = UkA reduces to UA/ng which increases with- 
out limit as the B-component density n g falls to zero for 
fixed UA, and in fact the continuous momentum distribu- 
tion [L /2'K)nkA reduces to a representation of the ideal 
Bose gas distribution NAS{k) as ng 0. UkB has the 
same behavior, with A and B interchanged. 

Two-particle density matrices, pairing, and ODLRO: 
The two-particle A-component and B-component density 
matrices P2Aa{xi,X2;x[,x'2) and /02Bs(yi, 2/2; yi, 2/2) can 
also be evaluated in closed form by generalization of ([3]) 
and the derivation for the one-component FTG gas in 
For component A 

P2Aa{xi,x2;x[,x'^) = Na{Na - 



X J ^'o{xi,X2,X3,- ■■ ,XNa;Y) 

X '^a{x[,x'2,X3,--- ,XNA;Y)dx3---dxNAdY (4) 

Using ■^o{xi,X2,X3,- ■ ■ ^xna'jY) = Ilf=^iSgn{xi - 
?/j)sgn(x2 — Uj) one finds in the thermodynamic limit 
P2AAixi,X2;x[,x'2) = n^e^"^ (^1-^2+23-24) where zi < 
Z2 < Z3 < Z4 are the arguments (xi ) in 

ascending order. If xi < X2 < 

P2AA = ? 



-2n_B — a;2 I p — 2nB |x J -Xj 



Xi < X2 then 
Generalizing the 
argument in [12] one sees that if the variable pairs 
(xi,X2) and {x[,x'2) are separated to arbitrary distance 
while keeping \xi — X2I and | I fixed, then p2AA 

remains constant, signalling AA-pair ODLRO associ- 
ated with a leading term Xiaauia{xi,X2)uia{x'i,X2) in 
the spectral representation of p2AA, with eigenfunction 
uia{xi, X2) = Ca e-^"^!^!-^^', normalization constant 
Ca ~ \/^ngjL, and macroscopic eigenvalue Xiaa = 
n\/C\ — TtA^ Al'^ng. There is a BEG-BCS crossover 
from AA-pair BEG when ng ^ ua and the range of 
uiA is <C y/nA implying tightly bound AA pairs, to AA- 
pair superconductivity when ng nA and the range 
of uiA is ^ ^/nA implying extended and strongly over- 
lapping AA Gooper pairs. Since /32BB(yi, 2/2; 2/2) ex- 
hibits the same behavior with A and B interchanged. 



one concludes that when ng ^ nA there is coexistence 
of BEG of AA pairs and superconductivity of BB pairs, 
and when ng <^ nA The A A and BB pairing is a purely 
off-diagonal phenomenon, both in the case of supercon- 
ductive ODLRO with weakly bound Gooper pairs and 
in the case of BEG of tightly-bound pairs. There is 
no diagonal AA, BB, or AB order; the pair distribu- 
tion functions DAAixi,X2) = n'^^ p2AAixi, X2; Xi, X2), 
Dgg{yi,y2) = ng'^p2gg{yi,y2;yi,y2), and DAB{x,y) = 
{nAng)^^ P2ABix, y] X, y) are all constant, as is most eas- 
ily seen by noting that (a) our system of A-bosons and 
B-bosons with FTG AB attraction maps to a mixture of 
ideal A-Bose and B-Bose gases with no AB interaction, 
and (b) diagonal density matrix elements are invariant 
under mapping via Eq. ([2]). This generalizes the previ- 
ous result for the pure FTG gas, where there is super- 
conductive ODLRO but the pair distribution function is 
constant 



The two-particle AB density matrix is also of interest. 
By a derivation paralleling that for p2AA and p2BB one 
finds in the thermodynamic limit 



P2AB{x,y;x' ,y') = nAngsgn{x - y)sgn{x' - y') 



X e 



-2nB\x~x' \ -2nA\y-y' \ 



(5) 



Suppose that x' = x + d and y' = y + d. Then 
P2AB — n^nBe-^^"^^"^'''' which vanishes exponentially 
as d ^ 00. It follows that there is no AB-pair ODLRO 
in spite of the strong AB attraction. It is informative in 
this connection to compare and contrast two cases (a) two 
ideal Bose gases A and B with no AB interactions, and 
(b) the present case, two ideal Bose gases A and B with 
FTG AB interactions. In case (a) there is complete BEG 
of both components A and B, the many-body ground 
state is a trivial constant, p2AA, P2BB, and p2AB are also 
constant, and hence all three of these density matrices 
exhibit ODLRO. However, the AA and BB ODLRO is 
a trivial consequence of the trivial ground state struc- 
ture and more generally, follows from the ODLRO of piA 
and piB [lij and implies no true AA or BB pairing cor- 
relation. In the present case (b), there is no BEG of 
component A or B and hence no ODLRO of piA or pis, 
but there is ODLRO of both p2AA and p2BB, implying 
both AA and BB pairing, since uia and uig have finite 
range. On the other hand, p2AB has no ODLRO. 

Quantum phase transition: Suppose now that in addi- 
tion to the odd-wave AB interaction of FTG form, there 
is also an even- wave AB interaction of LL delta function 
form [ill, Ve{xi - yj) = gAB^^Xi - yj) with cjab > 
(repulsive interaction). Recall that before passing to the 
infinitely narrow well limit xq — > 0, Vq — > 00 of the FTG 
interaction, the ground state in the case of no even-wave 
interaction has nodes at AB contact due to the inter- 
nal wave function sinK^Xi — yj). These nodes kill the 
even-wave interaction, so that the ground state in the 
presence of the even- wave interaction is the same as that 
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in the absence of an even-wave interaction, which has 
energy zero. Any state without AB contact nodes has 
positive energy in the presence of the even-wave repul- 
sion, since then the odd-wave FTG interaction is killed; 
the system develops AB contact nodes spontaneously in 
the presence of the FTG odd-wave attraction in order 
to lower its energy by killing the even- wave repulsion. 
Now suppose that there is no even-wave interaction, but 
the odd-wave attraction is weakened by carrying out the 
limit .To and Vq ^ cx) in such a way as to produce 
a large but finite negative odd-wave scattering length, 
— oo < QABo < 0. So far the FTG limit where kxq = f 
and ttABo — —oo has been assumed, but more generally 



2x0 



as 



if lABo is negative and finite, then kxq = ^ Tr\aAB I 
xq — > [1, . The ground state in that case maps 
via to a model ground state ^mo consisting of two 
Bose gases A and B with no AA or BB interaction but an 
even- wave AB interaction of LL form g'^gd{xi — ijj) with 
9'ab — I lA^ABo\ and positive energy. If one general- 
izes further by adding a nonzero even-wave interaction 
gAB^ixi — Dj) again, then so long as the ground state 
is retained unchanged with AB contact nodes, the even- 
wave interaction will be killed and will have no effect. 
However, if g^B < ^'ab ~ I ^AO'ABo\l then a ground 
state with no AB contact nodes has lower energy, since 
then it is the odd-wave interaction which is killed, instead 
of the even-wave interaction. It follows that there is a 
quantum phase transition between a phase with no AB 
contact nodes, which has lower energy when gAB < q'ab^ 
and a phase with AB contact nodes, which has lower 
energy when gAB > Qab- This is very similar to the 
ferromagnetic-antiferromagnetic phase transition in the 
ID spinor Fermi gas with both even and odd-wave inter- 
actions l20il2H. 

Stronger AB attraction: It was pointed out above that 
if — oo < OABo < then kxq — ^ ~ Ti\a^^g Xq — > 0, 
or without the absolute value signs kxq — ^ ' 



In 



analogy with the case of fermions discussed recently [24] , 
for an interaction of FTG form the same relation (with- 
out absolute value signs) holds if oabo > 0, in which 
case the model state 5'7\/o generated by the mapping 
(PI) has an attractive LL AB interaction gAB^ixi — yj) 
with gAB = —fi^/i^o,ABo < and negative energy. The 
N = 2 model ground state is bound, ^'oa/ = e~l^~^l/°-*^° 
with energy £^0 — ~^r~T — ^ and the corresponding ac- 

tual physical state is ^'o = sgn(a; — i/)e~l^^^^l/'^^^°, with 
the same energy. The exact solution for both Na > 2 
and Nb > 2 is not known, but if has Bijl-Jastrow 
form ^Mo = ^n.=l nfii fi\x^ - Vjl) where 91 is a nor- 
malization constant and f{\xi — yj\) ~ Q^\^i~Vi\l°-ABo fQj. 

~ yj\ ^^fnAJiB)'^ , then in the neighborhood of 
each Xi — yj — the wave function reduces to that of an 
AB dimer. The behavior of '^mo as Xi recedes to dis- 
tances Xi ^ (y'ri^ns)^^ from yj will be controlled by 
proximity of Xi to other B-particles, not yj. If /(^) X 



for ^ » {y/nAnB)~^ where A is some nonzero constant, 
then = A^^^ . The physical state ^0 generated by the 
mapping ^ then has reduced density matrices reducible 
to ID integrals as previously. However, in the thermody- 
namic limit P2AB reduces to the previous expression (O, 
with no AB-pair ODLRO. A better approximation to ^'o 
might reverse this conclusion, so this model deserves fur- 
ther study. 

Outlook: In view of the surprising properties of this 
model, experimental study of mixtures of weakly inter- 
acting A-bosons and B-bosons in tight waveguides with 
strong ID AB attractions 0] induced by a p-wave AB 
Feshbach resonance should be fruitful. In particular, 
one could look for the predicted quantum phase transi- 
tion. 

I thank Peter Reynolds and Gregory Astrakharchik for 
helpful comments on drafts of this work. 
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